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Abstract 

In dynamic atomic force microscopy, nanoscale properties are encoded in the higher harmonics. Nevertheless, when gentle interac- 
tions and minimal invasiveness are required, these harmonics are typically undetectable. Here, we propose to externally drive an 
arbitrary number of exact higher harmonics above the noise level. In this way, multiple contrast channels that are sensitive to 
compositional variations are made accessible. Numerical integration of the equation of motion shows that the external introduction 
of exact harmonic frequencies does not compromise the fundamental frequency. Thermal fluctuations are also considered within the 
detection bandwidth of interest and discussed in terms of higher-harmonic phase contrast in the presence and absence of an external 
excitation of higher harmonics. Higher harmonic phase shifts fijrther provide the means to directly decouple the true topography 
from that induced by compositional heterogeneity. 



Introduction 

It has long been recognized in the community that higher 
harmonics encode detailed information about the non-lineari- 
ties of the tip-sample interaction in dynamic atomic force 
microscopy (AFM) [1-5]. Physically, non-linearities relate to 
the chemical and mechanical composition [6] of the tip-sample 
system and imply that higher harmonics can be translated into 
conservative and dissipative [7] nanoscale and atomic prop- 
erties [8]. Furthermore, conventional dynamic AFM can already 



reach molecular [9,10], sub-molecular [11] and atomic [12,13] 
resolution in some systems. Thus, the simultaneous detection 
and interpretation of multiple higher harmonic signals while 
scanning [14] can lead to spectroscopy-like capabilities [15,16], 
such as chemical identification, with similar or higher resolu- 
tion [5,17,18]. The higher harmonic approach however, and 
particularly in other than highly damped environments [19,20], 
requires dealing with the recurrent challenge of detecting higher 
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harmonics [1,3:21,22]. Higher harmonics are a result of the 
non-linear tip-sample interaction in the sense that the inter- 
action effectively acts as the driving force of each harmonic 
component [7]. Accordingly, relatively high peak forces, of the 
order of 1-100 nN, are required [22,23] to excite higher 
harmonics above the noise level. In order to address this issue, 
in 2004 Rodriguez and Garcia [23] proposed to drive the second 
higher flexural mode of the cantilever with an external drive. In 
this way, and by driving with sufficiently small (sub- 
nanometer) second mode amplitudes, the first mode amplitude 
[24] or frequency [17] can be employed to track the sample in 
amplitude or frequency modulation (AM and FM), respectively. 
The second mode can then be left as an open loop for high 
sensitivity mapping of compositional variations [25] or as a 
closed loop, in which case the tip-sample stiffness ^ts can be 
computed [17,26]. More recently, the multitrequency AFM ap- 
proach has been extended to employ three flexural modes [27] 
and/or simultaneous torsional modes [28], for which, typically, 
the frequency and mode under consideration are externally 
excited [24]. In summary, FM and/or AM feedback systems can 
be employed in one [29], several [27] or all of the modes under 
consideration in order to quantify properties on the nanoscale 
through observables [30] while simultaneously enhancing sensi- 
tivity and throughput [31]. The dynamics in the multifirequency 
approach, however, might lead to extra complexities in the 
analysis, acquisition and interpretation of data [31,32]. For 
example, recent studies [31] show that multiple regimes of 
operation might follow depending on the relative kinetic 
energy between the higher mode of choice and the fimdamental 
eigenmode [31,33]. 

Here, exact multiple harmonics of the fundamental drive 
frequency are externally excited above the noise level to open 
multiple contrast channels that are sensitive to compositional 
variations. The focus is on amplitude modulation (AM) AFM, 
in which the fundamental amplitude A\ = A tracks the sample 
as usual. For standard cantilevers the eigenmodes are nonhar- 
monic [29]. That is, the natural resonant frequencies of the 
cantilevers are not integer multiples. Furthermore, these natural 
frequencies relate to the geometry and mechanical properties of 
the cantilever [34]. The practical implication is that it is only 
easy to induce large oscillations at the frequencies that coincide 
with these natural frequencies. Nevertheless the tip-sample 
coupling always occurs via harmonic frequencies. This is 
because a periodic motion always implies that there is a funda- 
mental frequency and that all other higher frequencies are 
integer multiples of the fundamental [35]. The implication is 
that externally introducing frequencies other than harmonic 
frequencies could induce a fundamental sub-harmonic 
frequency [24,35]. In short, the incommensurability between 
external drives in the standard multifrequency approach implies 



that the cantilever motion is not exactly periodic relative to the 
fimdamental drive and that a sub-harmonic excitation typically 
follows [32]. Furthermore, simplifications in eigenmode 
frequency shift theory [36] might lead to inconsistencies [37]. 
This issue becomes more prominent when dealing with third or 
higher eigenmodes [27,38], for which the theory is now 
emerging [31]. The introduction of exact harmonic external 
drives keeps the fundamental frequency intact and the analyt- 
ical expressions are simplified by orthogonality. Furthermore 
2{N-\) observables, i.e., higher harmonic amplitudes and 
phases, are made available even with peak forces no higher than 
200 pN, as they are required [25,39] for high resolution and 
minimally invasive imaging of soft matter. Thermal fluctua- 
tions are also considered here in order to establish a possible 
loss of contrast due to fundamental sources of noise. It is also 
shown that true topography and apparent topography, which is 
induced by chemical heterogeneity, can be decoupled at once by 
monitoring the phase contrast of higher harmonics. 

Results and Discussion 

Consider the equation of motion of the /nth eigenmode 

»(m) y(m)«(m) 

where Q(m)y and Z(„) are the spring constant, quality 
factor, natural frequency and position of the wth eigenmode. 
The term Fj) stands for the external driving force 

N 

^D(0=X^O«cos(na)r), (2) 

n=l 

where the subscript without brackets, n, indicates the harmonic 
number. Note that here m„ = nm, where co is the fundamental 
drive frequency set near mode m = 1, i.e., o) = (aj = a^^iy The 
term Fjs is the tip-sample force, which is a fimction of both the 
tip-sample distance, d, and velocity, d. Here however, we 
focus on conservative forces since these are present even with 
gentle interactions. Hence we can write i^tsW- Since the higher 
harmonic amplitudes here are externally excited, the number of 
harmonics A' that is to be monitored can, in principle, be arbi- 
trarily chosen up to the limits of frequency detection, i.e., of the 
order of MHz, without compromising detection. The main 
constraint is that the number of higher modes, M, that is to be 
considered needs to be consistent with the number of higher 
harmonics N that are to be analysed [22]. For simplicity, we 
consider M = 2 and N= 1 0 in the numerical analysis without 
loss of generality. For clarity we emphasize that M is the 
number of modes and N is the number of harmonics taken into 
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consideration in the analysis in this work. A particular mode or 
harmonic is referred to in lower case, i.e., morn respectively. 

The nth harmonic velocity z„ is 



If the «th drive Fo„ is zero, then 



(9) 



^n[t) = -A^n sin(ncD«-(|)„). 



(3) 



Multiplying Equation 1 by Equation 3 and integrating over a 
cycle results in 



nk 



+71- 



Q 



(4) 



where 1 is assumed when no subscripts are given. The relation- 
ships (co/co(„))^ = k/k(^„j and Q/Q(m) = co/a>^„^ [7] have been 
employed in Equation 4 and it has been assumed that the funda- 
mental drive frequency co is set near 05(1). Furthermore, in Equa- 
tion 4 ^(ot)„ and A„ are the amplitudes of the «th harmonic that 
correspond to the position of mode m, i.e., Z(„), and to the 
absolute position of the tip, i.e., z, respectively. Also 



B 



(m)n = COS <t)„ sin ^^^j„ - sill (|)„ COS (|)(^)„ , (5) 



^{m)n ~ COS(|)„ COS(t)(^^-)„ + sin(t)„sin(|)(^)„ , (6) 



u=t+T 



(7) 



Equation 9 is the energy transferred to the nth harmonic of the 
cantilever through the tip-sample interaction. It should be noted 
that this is consistent with a conservative tip-sample force 
Fis{d) since the energy is provided during each cycle by the 
external driving force(s). The quadratic dependence of the 
energy E„ on nA„ is of particular relevance for the detection of 
higher harmonics. First, Equation 9 implies that for a given 
amplitude A„ the transfer of energy E„ scales quadratically with 
the harmonic number. This explains why for sufficiently large 
n, higher harmonics are typically undetectable. Second, the 
proportionality between E„ and in Equation 9 explains why 
for higher harmonic amplitudes to be detected, the interaction in 
Equation 7 needs to be considerably large, even when n is not 
necessarily very large. 

From Equation 8 it follows that A„ can be set to any arbitrary 

value by increasing Fq„, even if there is no tip-sample energy 
transfer, i.e., E„ = 0. The higher harmonics for the free 
cantilever are termed Aq„. This case corresponds to a free 
cantilever oscillating sufficiently high above the sample 
(A/Aq = 1) as illustrated in Figure 1 (circles). The data has been 
acquired by numerically solving the simultaneous equations in 
Equation 1 for the first two flexural modes, i.e., M= 2, and for 
7V= 10. Furthermore, since only long range attractive forces are 
of interest here, the tip-sample force is simply [23] 



RH 



Fts {d) = y with d> Uq , 



6d 



(10) 



where (])(„■)„ and (|)„ are the phase shifts of the nth harmonic that 
correspond to the mth mode position and the absolute position, 
z, respectively, and E„ is the energy involved with the nth 
harmonic tip-sample interaction. Near the modal frequency 
tt'(m) only the with mode significantly contributes to the inter- 
action and B(„)„ 0 and C(„)„ 1 in Equation 4. This approxi- 
mation has been currently employed in the literature [6]. Never- 
theless, far from the modes, these terms might not be zero. To 
allow for simple analytical formulae and ease the qualitative 
interpretation we consider the harmonics close to the modes 
only [6]. Then 



(8) 



where R is the tip radius. His the Hamaker constant and is an 
intermolecular distance (qq = 0.165 nm throughout and in all the 
data here, we consider d> oq throughout). It is relevant to note 
that the Hamaker constant depends on the tip and sample in the 
sense that its value is determined by the atomic composition or 
chemical elements that compose the tip and the sample [40,41]. 
For this reason, in this work we will employ the terms chem- 
istry, Hamaker and tip-sample composition or chemistry 
interchangeably. The common parameters in this work are 
k = 2 N/m, Q = 100, co = 2;f70 kHz and i? = 7 nm, i.e., they 
correspond to commercially available standard probes for AM 
AFM. Furthermore, in Figure \, H = 6.2 x 10~19 J, i.e., it is 
close to that calculated for materials such as polystyrene or 
fused quartz [40]. The parameters for the second mode have 
been obtained with the above formulae [7]. The modal Irequen- 
cies 1 and 2 are shown with dashed lines. The phase shifts 
are shown in the vertical axis in Figure 1 for each harmonic. 
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Figure 1 : Phase shifts t!f„ of higher harmonics, including the fundamental shift <t)i, when W = 10 external harmonic drives are introduced. The values 
(jip are shown for a free oscillating cantilever (circles). For the free cantilever the separation is » AIAq - 1 . Then the cantilever is gently interacting 
(peak forces smaller than 20 OpN) with the surface, i.e. , AIAq < 1 , while the free higher harmonic amplitudes Aq„ are set to 1 (squares) and 1 00 (trian- 
gles) pm. 



The actual harmonic amplitudes that resulted when inter- 
acting are not shown, instead A„ ~ Aq„ is given throughout. The 
case of a free cantilever (circles) shows that the fundamental 
phase shift is exactly 90 degrees as expected while the higher 
harmonic phase shifts (« > 1) lie either close to 180° or to 
0°. This is in agreement with Equation 8 when E„ ~ 0 since then 

where the approximation Fq,, ~ k^^yi^AQn (near m) has been 
employed. Also from Figure 1 (circles) it follows that for a free 
cantilever, and when n is higher than the modal frequency 
(close to a given mode and for n > 1), ~ 180°. When n is 
lower than the modal frequency (})„ = 0°. This is true irrespec- 
tive of the value of^o«- When the tip is allowed to interact with 



the sample E„i^O and, from Equation 8, the phase shift is 
affected by the interaction. Nevertheless, the weight of the 
driving force, i.e., the first term in Equation 8, increases with 
increasing Fo„, or Ao„, and then the sensitivity of (()„ to 
might be compromised. This is confirmed in Figure 1 by 
allowing a gentle interaction, i.e., ^oi = Aq = 4 nm and 
A/Af) ~ 0.9 (also Figure 2 and Figure 3), and monitoring 
(|)„when Ao„ = 1 pm (squares) and Aq,, = 100 pm (triangles). 
Wheny^o,, = 100 pm (triangles) all ([)„ remain close to 180° or 
0°. A shift in phase, i.e., from 180° to 0°, is observed for « = 2 
only. While these jumps of nearly 180° might be of interest they 
are ignored from now on. The reader can refer to recent works 
that discuss multiple regimes of operation in bimodal AFM 
[31,33]. It follows that variations in Hamaker are not detected 
by higher hamionic frequencies when Aq„ = 100 pm. When Aq„ 
= 1 pm (squares), however, the values of (()„ are not exactly 
180° or 0° for some n. Thus, the values <()„ are now sensitive to 
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Figure 2: Phase shift analysis, in which the contrast in the higher harmonic phase shifts Aijj^ = abs(ij>„(H2) - i|)n(Hi)), n = 2-10, which is induced by 
variations in the Hamaker constant H is shown. The variation in /-/ is H2 - Hi = 1 .0 x 10"''^ J, where H2 = 1 .2 « 10"^^ J, and effectively corresponds to 
variations in chemistry only. Results are shown when higher harmonic amplitudes Aq,, of 1 (circles), 10 (squares) and 100 (triangles) pm are intro- 
duced. Peak forces are smaller than 200 pN throughout. 
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Figure 3: Phase shift analysis, in which the contrast in the higher harmonic phase shifts Acji^ = abs( ^„ (H2) - it„ (Hi)) for n = 2-1 0 results only from 
variations in the Hamaker constant, H, or In the chemistry. The variations of Hare H2 - Hi = 0.2 x 10"''^ J for H2 = 0.4 » 10"''^ J (circles), 0.8x1 0'''^ J 
(squares) and 1 .2 x 1 0"''^ J (triangles). These variations Induce variations In peal< force of 29 (circles), 8 (squares) and 3 (trinagles) pN. 



the Hamaker values or tip-sample forces. The peak forces were 
140 pN (circles) and 160 pN (triangles) respectively. 

The loss of phase sensitivity to Hamaker variations with 
increasing Aq„ is flirther corroborated with the use of Figure 2 
and by varying the Hamaker values from //i = 0.2 x 10~'^ J to 
H2 = lA X 10~" J, and setting = 1 P™ (circles), ^on =10 P™ 
(squares) and Af)„ =100 pm (triangles). This range of H 
is characteristic of materials interacting in ambient conditions 
[40]. The >'-axis stands for the contrast in higher harmonic 
phase A(()„ = abs((t)„(//2) ^ 't'«(^l))- We consider that varia- 
tions, for which A(|)„ > 0.2° lie above the noise of the instru- 
ment and can potentially be detected. The corresponding varia- 
tions in peak forces were 63, 47 and 79 pN respectively. The 
sensitivity of A^^ is clearly controlled by the chosen values of 
Aq„ For example, if Aq,, = 100 pm then A(()„ < 0.2° throughout. 
If ^On = 1 or 10 pm, however, then A(()„ > 0.2° at least for some 
n. In particular, if Aq„ = 1 pm then A(|)„ > 0.2° for all n. This 
implies that all the externally excited higher harmonics act as 
simultaneous contrast channels that are sensitive to Hamaker, or 
chemical, variations. 

In Figure 3 the sensitivity of A(()„ when Ao„ = 1 pm is tested by 
varying H (a) from Hi = 0.2 x IQ-I^ J to //2 = 0.4 x IQ"!'^ J 
(peak force variation of 29 pN, circles), (b) from 
Hi = 0.6 X 10"'^ J to //2 = 0.8 X 10"15 J (peak force variation 
of 8 pN, squares) and (c) from //[ = 1.2 x 10"'^ J to 
7/2 = 1 -4 X 10~'^ J (peak force variation of 3 pN, triangles). The 
shifts A(()„ are larger than 0.2° for all n provided the variations 
in peak force are large enough (circles). If the variations in the 
peak force are sufficiently small then A(|)„ > 0.2° for some n 
only. Also, it can be deduced by inspection that, in general, A^^ 
escalates with variations in peak force and changes non-linearly 
with variations in Hamaker since Hj - Hi = 0.2 x 10"'^ J 



throughout in the figure. In fact, from Figure 3, the total contri- 
butions to the phase shift calculated as the sums SA(|)„ (n = 1-9) 
are 1 19.8, 19.3 and 5.4° and decrease with decreasing the varia- 
tions in peak force, i.e., 29, 8 and 3 pN, respectively. 

It is also interesting to note that the source of variations in peak 
force with variations in Hamaker H (Equation 10), i.e., van der 
Waals forces, relates to variations in the distance of minimum 
approach, d^, with variations in H. To be more specific, d^, 
increases with increasing H. For example, in the simulations, by 
varying H from Hi = 0.2 x IQ"!'' J to //2 = 1-4 x lO"'** J the 
variation is Ad^ ~ 0.83 nm. This would experimentally result in 
a chemistry-induced apparent topography of approximately 
AZ(. ~ 0.83 nm. In standard AM AFM, in which a single 
frequency is externally excited, this apparent topography cannot 
be distinguished from true topography in the presence of 
conservative forces only (Figure 4). A true topography can only 
be reconstructed from AM AFM results, if there is a variation in 
topography only (Figure 4a). This means that the composition 
of the sample is homogeneous throughout. In particular, the 
above discussion indicates that variations in H, or chemistry 
alone, produce variations in apparent topography in AM AFM, 
for which Az^ > 0 nm (Figure 4b). The excitation of higher 
harmonics, however, provides experimental observables 
to differentiate between the two cases. Namely, the true 
reconstructed topography results only if A(|)„ = 0° for all n. That 
is, if A<|)„ > 0°, even for a single n, there is a contribution to 
apparent topography induced by chemistry or other composi- 
tional variations. 

Thermal noise and higher harmonic external 
drives 

As stated in the introduction, it has long been known that under 
ambient conditions higher harmonic amplitudes might be too 



272 



Beilstein J. Nanotechnol. 2014, 5, 268-277. 



TRUE TOPOGRAPHY 




Figure 4: (a-c) Illustration of a cantilever oscillating above a surface 
and recovering tine true fieigint Az^ = h when there are no composi- 
tional heterogeneity or chemical variations, (d-f) Topographical varia- 
tions AZc > 0 nm induced by chemical or another compositional hetero- 
geneity. The two cases can be decoupled by noting that it is only a 
compositional heterogeneity, if the phase shifts of higher harmonics, 
A(|)n, are non-zero. 

small to be detected [3,15,42]. This is particularly true when 
monitoring higher harmonics and simultaneously applying 
gentle tip-sample forces [23]. In liquid environments, however, 
the second harmonic amplitude might be large enough [43] to 
be recorded to map the properties even of living cells [44,45]. 



Still, even in highly damped environments, harmonic ampli- 
tudes rapidly decrease with increasing haiTnonic number paitic- 
ularly when imaging soft matter [6,15,46]. The main discussion 
above has focused on externally driving higher harmonics to 
amplitudes that could be experimentally detected. Then, once 
these amplitudes are sufficiently high, the phase shifts A(|)„ have 
been employed to map the composition through variations in 
the tip-sample Hamaker constant, H, in Equation 10. In this 
section, the presence of thermal noise is discussed with respect 
to the contrast in amplitude AA„ and phase A(|)„ in the presence 
and absence of external drive forces at the higher harmonics 
frequencies. 

First an example of the magnitude of the harmonic amplitudes 
and respective phase shifts that would result when higher 
harmonics are not externally excited is given (Table 1). In order 
to sense long-range forces only, the cantilever is driven with 
relatively small amphtudes, i.e., Aq = 4 nm and A/Aq = 0.9 as in 
the examples above. The harmonic amplitudes A„ are given in 
pm. Two examples for the amplitude response are shown, one 
for amplitudes resulting from Hi = 0.2 x 10~" J (top row) and 
one for //2= 1.4 x 10~'^ J (second row). For //i, ^2 is approx. 4 
pm whereas A^ and A^ are approx. 1 pm. All other higher 
harmonics lie below 1 pm. For H2, Aj is approx. 3 pm and all 
other higher harmonics have values below 1 pm. The differ- 
ence in amplitudes KA„ = A„{H2) -A„(Hi) that results from the 
variation in H is also given in the table. Only the second 
harmonic results in variations above 1 pm. Practically, these 
results imply that while higher harmonic amplitudes depend on 
the value of the Hamaker constant, or sample composition, the 
amplitude values are typically in the order of 1 pm or fractions 
of a pm. This is also true for variations in higher harmonic 
amplitudes AA„. The corresponding phase shifts (])„ and varia- 



Table 1 : Harmonic amplitudes and the corresponding phase shifts ^„ that result from Hamaker values of Hi = 0.2 x 1 0~^^ J and 
H2 = 1 .4 X 10"^^ J. The differences in amplitudes AA„ and phases are also shown. A single external drive force has been employed (the funda- 
mental frequency) and no thermal noise has been allowed. 
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tions in phase shifts A(|)„ are also shown in Table 1 for Hi and 
H2. These are of the order of a hundredth of a degree or less 
except for sufficiently high harmonic numbers, i.e., n = 9 and 
10. The amplitudes for these higher harmonics, however, are of 
the order of tens of femtometers or less. 



J ^ 

f„-Afl2 kU [l-CO^ /co2„)]' +[co/(cO(„)e(„))]' 



1 



df. 



(15) 



Thermal fluctuations are a fundamental source of intrinsic noise 
in atomic force microscopy [47]. Thus, while other sources of 
intrinsic and extrinsic noise should be acknowledged and might 
be present in a given experiment, thermal fluctuations are 
analyzed next in terms of their effects on amplitude and phase 
shifts. This should provide a measure of the impact of thermal 
noise on the enhanced contrast reported in this work (Figure 2 
and Figure 3). Other technical issues such as tilt and probe 
geometry have also been ignored for simplicity since these typi- 
cally involve a correction factor [48]. As in the work of Butt 
and Jaschke [47], the equipartition theorem is employed to esti- 
mate the thermal noise present in a given mode. However, since 
higher harmonics are discussed here, particular emphasis should 
be given to the noise at the frequencies of interest, i.e., at exact 
harmonic frequencies, and the noise in the detection band- 
widths of interest. Then, the thermal noise power APxn(4/) in 
the detection bandwidth of interest. A/, can be defined as 

/„+A//2 

APtn(A/)-2 J Gtn|/^zf| df, (12) 
/„-A//2 

where TN stands for thermal noise, f„ is the frequency of 
interest (cy„ = 2nf„), that is the frequency of a particular 
harmonic n, G™ is the power spectral density due to thermal 
noise, and |/^zfP is the modulus of the squared transfer func- 
tion of a particular mode m of position z„, relative to thermal 
force Fj^. If Gjn is assumed to be constant for the bandwidth 
of interest in AFM experiments, i.e., f = 10^-10^, it follows 
from Equation 1 that the thermal energy in a given mode m, by 
invoking the equipartition theorem, is 



rknT 



^(»i)'^TN j 



f ''^•(13) 



where here T = 300 K throughout, /(„) is the natural resonant 
frequency of mode m in Hz and df= (/(„)/a)(„))fifcu. Then 



(14) 



From Equation 13 and Equation 14, the thermal noise power in 
the detection bandwidth of interest, APxn(A/)> is found to be 



Finally, the associated amplitude due to thermal noise ^xN iii 
the detection bandwidth A/ is 



^tn(A/) = V2APtn(A/). 



(16) 



It should be noted that ^xN gives the contribution of thermal 
noise to the amplitude of a given mode m only. Each modal 
contribution of thermal noise to the amplitude should be calcu- 
lated separately for each frequency in the formalism developed 
here. A driving force, -Fxn, can also be associated to thermal 
noise and the respective amplitude, ^xN> (Equation 16) through 
a standard expression [49] 



CO 
«>(m) 



CO 



(17) 



Equation 17 gives the effective drive force Fym due to thermal 
fluctuations that should be expected for a given detection band- 
width A/ and a given mode m. Since the upper boundaries for 
noise will be considered here, the phase of the thermal noise 
signal has been set to be in quadrature with respect to the 
external drive, i.e., either the fundamental external drive or the 
higher harmonic external drives when these are present. Focus 
is now placed on the harmonics n = 1, 2, 3 (close to the funda- 
mental frequency of mode 1) and 6 (close to the fundamental 
frequency of mode 2), since these are sufficiently close to a 
given mode that only the contribution of thermal noise to the 
amplitude from a single mode needs to be considered. This 
simplifies the following discussion. 

In Table 2 the amplitudes ^xNn forces Fj^„ calculated for 
three different values of detection bandwidth A/ (5 kHz, 2 kHz 
and 0.2 kHz) are shown for « = 1, 2, 3 and 6. The values have 
been computed with the use of Equation 16 and Equation 17, 
with frequencies centered at the harmonic frequencies y^, for a 
given detection bandwidth A/ It is interesting to note that^xNl 
lies between 44 and 19 pm for the three choices of detection 
bandwidth. These values are in agreement with those expected 
from an analysis that implies that all the thermal noise is 
centered exactly at resonance [47]. This is because the Q factors 
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are relatively high (Qi = 100 and Qj = 600). The values of the 
thermal-noise amplitude expected at harmonics 2, 3 and 6 
however are of the order of 0.1-1.0 pm. 



The effects that the thermal noise amplitudes in Table 2 have on 
the enhanced contrast reported in this work have been analyzed 
by adding the associated thermal noise forces, also shown in 
Table 2, to the equation of motion in Equation 1 . The discus- 
sion below focuses on the values obtained for A/= 2 kHz in 
Table 2 since this is a detection bandwidth of practical rele- 
vance in standard AFM experiments [50]. 

The sensitivity of the phase shift to noise and signal can 
be defined here, and for the purpose of phase shifts in 
AM AFM, as follows. First assume that noise is allowed 
according to Table 2 (A/= 2 kHz) for a given value of the 
Hamaker constant, H. Here both Hi = 1.4 x 10"'^ J and 
Hj = 1.4 X 10~19 J have been used in the simulations. 
According to this, thermal noise alone should lead to a differ- 
ence in phase shift A (])„(//) = (()„(.4tn > 0) - (|)„(.4xn = 0) for a 
given value of H since there is an effective driving force i^xNn 
due to thermal fluctuations (Table 2). The average of for 
the two Hamaker values can be taken as the noise in the phase 
signal as follows 

A^„(TN)= ^"^ ^ ^"^ ^ ^, (18) 



PR(4>„)= J" (20) 
A(t)„(TN) ^ ^ 

Large values of PR result in a high sensitivity of the phase shift 
to the signal, whereas low values of PR indicate a sensitivity of 
the phase shift to noise only. Three cases are discussed, which, 
for simplicity, focus on harmonics 2, 3 and 6 only and on 
= 0, 1 and 10 pm. 

Case 1: First, no higher harmonic external drives are allowed, 
which implies that Aq„ = 0 in Equation 2 for n> \. This is the 
standard operational mode in dynamic AFM, in which a single 
external drive is employed. In this case we have PR = 0 
throughout (Table 3). 

Case 2: Higher harmonic external drives are allowed. In par- 
ticular, Aq„= I pm in Equation 2 for n> 1. This is the proposed 
mode of operation in this work. In this case we have PR > 1 
throughout but the exact value depends on harmonic number 
(Table 3). 

Case 3: Higher harmonic external drives are allowed. In par- 
ticular, Aq„ =10 pm in Equation 2 for n > 1. This is the 
proposed mode of operation in this work. When compared to 
case 2, however, the magnitudes of the external drives have 
been increased. In this case we also have PR > 1 throughout 
(Table 3). 



Table 3: The phase ratio for a given harmonic phase shift n, PR{^„), 
as defined by Equation 20 when 1) no higher harmonic external drives 
are allowed (Ao„ = 0) and when external drives lead to 2) Aq„ = 1 pm 
and 3)/Aon= 10 pm. 

PR(<t)2) PR(<t>3) PR(<|)6) 



case1:/Aon=0 0.00 0.00 0.00 

case2:/\on=1 pm 1.90 22.09 7.29 

case3:/\on=10pm 5.20 2.01 195.85 



Table 2: Amplitudes At^„ resulting from thermal noise for a7 = 1 , 2, 3 
and 6 and respective drive forces Ftno for detection bandwidths Af of 













"B. " 








Af 


'^TN1 




'^TN2 












[kHz] 


[pm] 


[pN] 


[pm] 


[pN] 


[pm] 


[pN] 


[pm] 


[pN] 


5 


62.23 


1.27 


0.42 


2.83 


0.17 


2.83 


0.42 


2.83 


2 


56.57 


1.13 


0.28 


1.70 


0.14 


1.70 


0.28 


1.70 


0.2 


26.87 


0.57 


0.08 


0.57 


0.04 


0.57 


0.07 


0.57 



where TN stands for thermal noise as usual and A(j)„(TN) stands 
for the difference in phase shift at harmonic n that induced by 
thermal noise alone. Next the signal is defined as the phase shift 
induced by variations in Hamaker alone 

Acp„ =(p„(^^2)-9«(^^l)- (19) 

Finally, a parameter that quantifies the sensitivity of the phase 
shift to noise and signal, the phase ratio PR((()„) can be defined 
from the ratio between Equation 19 and Equation 18: 



When looking at Table 3, one should recall that these are the 
upper-boundary values for noise since the phase of the thermal 
noise drives was set to be in quadrature. In summary. Table 3 
shows that the phase ratio PR increases when external drives are 
applied at a given exact harmonic fi^equency, i.e., when At)„ > 0. 
This is consistent with standard multifrequency operation, for 
which impressive results have already been obtained by exciting 
fi^equencies close to the resonant frequency of the second flex- 
ural mode [17,25,26]. In standard monomodal dynamic AFM, 
in which a single external drive is employed, the higher 
harmonics are excited by the tip-sample interaction according 
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to Equation 9. That is, energy needs to flow into tlie higlier 
harmonic frequencies in order to increase the amplitude signal. 
It is reasonable to assume that the increase in the sensitivity of 
the phase shift to the signal, i.e., the force, when external drives 
are applied is a consequence of energy both entering and 
leaving the given harmonic frequency of choice. That is, the 
fact that energy is supplied by the external drive at a given 
harmonic n implies that both positive and negative energy 
transfer might also occur at that frequency. Furthermore, when 
external drives are employed, this transfer occurs for a given 
phase shift that is now measured relative to the angle of the 
driving force. This is in agreement with the presence of the 
phase shift in Equation 8 and the absence of the phase shift in 
Equation 9 and might be related to the increase in the sensi- 
tivity of the phase shift to the tip-sample force as predicted 
here. 

Conclusion 

In summary, we have introduced a method that makes readily 
accessible an arbitrary number of exact higher harmonics by 
externally driving them with amplitudes above the noise level. 
Driving with exact higher harmonics does not introduce sub- 
harmonic fi^equencies to the motion and the amplitudes do not 
significantly decay when the interaction is gentle. Once higher 
harmonic amplitudes are accessible, one can also detect varia- 
tions in higher harmonic phase shifts. In this work, variations in 
sample composition, or chemistry, here modelled through the 
Hamaker constant, have been shown to lead to variations in 
higher harmonic phase shifts and amplitudes. In particular, vari- 
ations in the Hamaker constant of the order of 10^" J can in- 
duce higher harmonic phase shifts in the order of 10°. This is 
provided the higher harmonic amplitudes are small enough, i.e., 
about 1-10 pm. These small variations in phase shift would 
suffice to distinguish between metals such as gold, silver or 
copper [40]. Higher harmonic phase shifts also provide the 
means to decouple the true topography from an apparent topog- 
raphy, which is induced by compositional variations. Further- 
more this outcome should still be valid in standard bimodal 
imaging. Overall, the proposed approach, and variations, might 
ultimately fulfil the promise of rapid chemical identification 
with multiple contrast channels while simultaneously exerting 
only gentle forces on samples. Still it has to be acknowledged 
that, experimentally, it is expected that technical issues might 
arise from the multiple excitation of exact frequencies and from 
the set-up required to detect variations in higher harmonic 
phase. In particular, the set-up would require the generation of 
exact harmonic external drives to bring the harmonic ampli- 
tudes above the noise level while keeping them small enough to 
provide enough phase contrast. This last point is relevant since 
it has been shown that higher harmonic amplitudes should 
remain in the sub-lOO-pm range for the higher harmonic phase 



shifts to be significantly large, i.e., above 0.2°, in response to 
variations in the tip-sample force. On the other hand, an 
analysis of thermal fluctuation that exploits the equipartition 
theorem has also indicated that thermal noise should be of the 
order of 0. 1-1 .0 pm close to the higher harmonics modes. The 
implication is that the working amplitudes should lie in the 
range of 1 to 100 pm. The noise analysis has also shown that 
there is an increase in sensitivity of the phase shift to the 
tip-sample force when frequencies are externally excited. 
Nevertheless, ultimately, only experimental practice, implemen- 
tation, ingenuity and further theoretical advances in the field are 
to establish what the limits of this approach are. 

Acknowledgements 

This project was funded by Ministerio de Economia y Competi- 
tividad (MAT2012-38319). The artistic figure was designed by 
Maritsa Kissamitaki. 

References 

1. stark, R. W.; Heckl, W. M. Surf. Sci. 2000, 457, 219-228. 
doi: 1 0. 1 0 1 6/30039-6028(00)00378-2 

2. Stark, M.; Stark, R. W.; Heckl, W. M.; Guckenberger, R. 
Proc. Natl. Acad. Sci. U. S. A. 2002, 99, 8473-8478. 
doi:10.1073/pnas.122040599 

3. Sahin, 0.; Quate, C; Solgaard, O.; Atalar, A. Phys. Rev. B 2004, 69, 
16541 6. doi:1 0.1 103/PhysRevB.69.1 6541 6 

4. Diirig, U. New J. Phys. 2000, 2, 1-12. doi:10.1088/1367-2630/2/1/005 

5. Sahin, 0.; Quate, C; Solgaard, O.; Glesslbl, F. J. Higher-Harmonic 
Force Detection in Dynamic Force Microscopy. In Springer Handbook 
of Nanotechnology; Bhushan, B., Ed.; Springer Verlag: Berlin, 
Germany, 2010; pp 711-729. doi:10.1007/978-3-642-02525-9_25 

6. Payam, A. F.; Ramos, J. R.; Garcia, R. ACS Nano 2012, 6, 
4663-4670. doi:10.1021/nn2048558 

7. Garcia, R.; Herruzo, E. T. Nat. Nanotechnol. 2012, 7, 217-226. 
doi;10.1038/nnano.2012.38 

8. Garcia, R.; Proksch, R. Eur. Polym. J. 2013, 49, 1897-1906. 
doi:10.1016/j.eurpolymj.2013.03.037 

9. Fukuma, T.; Kobayashi, K.; Matsushige, K.; Yamada, H. 

Appl. Phys. Lett. 2005, 86, 193108-193110. doi;10.1063/1. 1925780 

10. Gotsmann, B.; Schmidt, C. F.; Seidel, C; Fuchs, H. Eur. Phys. J. B 
1998, 4, 267-268. doi:10.1007/s100510050378 

11. Gross, L.; IMohn, F.; Moll, N.; Liljeroth, P.; Meyer, G. Sc/ence 2009, 
325, 1110-1114. doi:10.1126/science.1176210 

12. Wastl, D. S.; Weymouth, A. J.; Glesslbl, F. J. Phys. Rev. B 2013, 87, 
24541 5-245424. doiil 0.1 1 03/PhysRevB.87.24541 5 

13. Glesslbl, F. J. Science 1995, 267, 68-71. 
doi; 1 0.1 1 26/science.267.51 94.68 

14. Kawai, S.; Hafizovic, S.; Glatzel, T.; Baratoff, A.; Meyer, E. 
Phys. Rev. 6 2012, 85, 165426-165431. 
doi:10.1103/PhysRevB.85.165426 

15. Stark, R. W. Nanotechnology 2004, 15, 347-351. 
doi; 1 0.1 088/0957-4484/1 5/3/020 

16. Nutter, C; Platz, D.; Tholen, E. A.; Hansson, T. H.; Haviland, D. B. 
Phys. Rev. Lett. 2010, 104, 050801-050804. 

dol:1 0.1 1 03/PhysRevLett.1 04.050801 



276 



Beilstein J. Nanotechnol. 2014, 5, 268-277. 



17. Kawai, S.; Glatzel, T.; Koch, S.; Such, B.; Baratoff, A.; Meyer, E. 
Phys. Rev. Lett. 2009, 103, 220801-220804. 
doi:10.1103/PhysRevLett.103.220801 

18. Hembacher, S.; GiessibI, F. J.; Mannhart, J. Science 2004, 305, 
380-383. doi:1 0. 1 1 26/science. 1 099730 

19. Basak, S.; Raman, A. Appl. Phys. Lett. 2007, 91, 064107-064109. 
doi:10.1063/1 .2760175 

20. Xu, X.; Melcher, J.; Basak, S.; Reifenberger, R.; Raman, A. 
Phys. Rev. Lett 2009, 102, 060801-060804. 

doi: 1 0.1 1 03/PhysRevLett. 1 02.060801 
21.Sahin, O.; Magonov, S.; Su, C; Quale, C. F.; Solgaard, O. 

Nat. NanotechnoL 2007, 2, 507-514. doi:10.1038/nnano.2007.226 

22. Gadelrab, K.; Santos, S.; Font, J.; Chiesa, M. Nanoscale 20M, 5, 
10776-10793.doi:10.1 039/c3nr03338d 

23. Rodriguez, T. R.; Garcia, R. AppL Phys. Lett. 2004, 84, 449-551. 
doi:10.1063/1. 1642273 

24. Lozano, J. R.; Garcia, R. Phys. Rev. Lett. 2008, 100, 076102-076105. 
doi: 1 0.1 1 03/PhysRevLett. 1 00.0761 02 

25. Patil, S.; Martinez, N. F.; Lozano, J. R.; Garcia, R. J. Mol. Recognit 
2007, 20, 516-523. doi:10.1002/jmr.848 

26. Martinez-Martin, D.; Herruzo, E. T.; Dietz, C; Gomez-Herrero, J.; 
Garcia, R. Phys. Rev. Lett. 2011, 106, 198101-198104. 

doi: 1 0.1 1 03/PhysRevLett. 1 06.1 981 01 
27.Solares, S. D.; Chawla, G. Meas. ScL Technol. 2010, 21, 125502. 
doi:1 0. 1 088/0957-0233/21 /1 2/1 25502 

28. Kawai, S.; Glatzel, T.; Koch, S.; Such, B.; Baratoff, A.; Meyer, E. 
Phys. Rev. B 2010, 81, 085420-085426. 
doi:10.1103/PhysRevB.81. 085420 

29. Proksch, R. AppL Phys. Lett. 2006, 89, 1 13121-1 13123. 
doi:10.1063/1 .2345593 

30. Herruzo, E. T.; Perrino, A. P.; Garcia, R. Nat. Commun. 2014, 5, 
No. 3126. doi:10.1038/ncomms4126 

31 . Kiracofe, D.; Raman, A.; Yablon, D. Beilstein J. Nanotechnol. 2013, 4, 
385-393. doi:10.3762/bjnano.4.45 

32. Stark, R. W. AppL Phys. Lett. 2009, 94, 063109-0631 1 1 . 
doi:10.1063/1 .3080209 

33. Chakraborty, I.; Yablon, D. G. Nanotechnology 20^3, 24, 475706. 
doi:10.1088/0957-4484/24/47/475706 

34. Steidel, R. An introduction to mechanical vibrations, 3rd ed.; John 
Wiley & Sons: New York, NY, USA, 1999. 

35. Tolstov, G. P.; Silverman, R. A. Fourier Series; Dover Publication: New 
York, NY, USA, 1976. 

36. Herruzo, E. T.; Garcia, R. Beilstein J. NanotechnoL 2012, 3, 198-206. 
doi:10.3762/bjnano.3.22 

37. Aksoy, M. D.; Atalar, A. Phys. Rev. B 2011, 83, 075416-075421. 
doi:10.1 103/PhysRevB.83.075416 

38.Solares, S. D.; Chawla, G. J. Appl. Phys. 2010, 108, 054901. 
doi:10.1063/1 .3475644 

39. Engel, A.; Mijller, D. J. Wat. Stnict MoL Biol. 2000, 7, 715-718. 
doi:10.1038/78929 

40. Israelachvili, J. Intermolecular and Surface Forces: Academic Press: 
Burlington, MA, USA, 1991. 

41 . Hamaker, H. C. Physica (Amsterdam) 1937, 4, 1058-1072. 
doi: 1 0. 1 0 1 6/S0031 -89 1 4(37)80203-7 

42. Stark, R. W.; Heckl, W. M. Rev. ScL Instrum. 2003, 74, 5111-5114. 
doi:10.1063/1. 1626008 

43. Preiner, J.; Tang, J.; Pastushenko, V.; Hinterdorfer, P. Phys. Rev. Lett. 
2007, 99, 046102-046105. doi:10.1 103/PhysRevLett.99.046102 



44. Dulebo, A.; Preiner, J.; Kienberger, F.; Kada, G.; RankI, C; 
Chtcheglova, L.; Lamprecht, C; Kaftan, D.; Hinterdorfer, P. 
Ultramicroscopy 2009, 109, 1056-1060. 
doi:10.1016/j.ultramic.2009.03.020 

45. Turner, R. D.; Kirkham, J.; Devine, D.; Thomson, N. H. 
AppL Phys. Lett. 2009, 94, 043901. doi:10.1063/1.3073825 

46. Raman, A.; Trigueros, S.; Cartagena, A.; Stevenson, A. P. Z.; 
Susilo, M.; Nauman, E.; Contera, S. A. Nat Nanotechnol. 2011, 6, 
809-814. doi:10.1038/nnano.201 1.186 

47. Butt, H.-J.; Jaschke, M. Nanotechnology ^995, 6, 1-7. 
doi:10.1088/0957-4484/6/1/001 

48. Helm, L.-O.; KappI, M.; Butt, H.-J. Langmuir 2004, 20, 2760-2764. 
doi:10.1021/la036128m 

49. Hu, S.; Raman, A. Nanotechnology 2008, 19, 375704-375714. 
doi: 1 0.1 088/0957-4484/1 9/37/375704 

50. Ramos, J. R. Desarrollos y applicaciones de la microscopi'a de fuerzas 
para el estudio de protei'nas y de celulas cancerosas. Ph.D. Thesis, 
Universidad Autonoma de Madrid, Madrid, 2013. 



License and Terms 

This is an Open Access article under the terms of the 
Creative Commons Attribution License 
( http://creativecommons.Org/licenses/by/2.0) , which 

permits unrestricted use, distribution, and reproduction in 
any medium, provided the original work is properly cited. 

The license is subject to the Beilstein Journal of 
Nanotechnology terms and conditions: 
( http://www.beilstein-joumals.org/bjnano ) 

The definitive version of this article is the electronic one 

which can be found at: 

doi:10.3762/bjnano.5.29 



277 



